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ABSTRACT 

Context. The source 4U 2206+54 is one of the most enigmatic high-mass X-ray binaries. In spite of intensive searches, X-ray pulsa- 
tions have not been detected in the time range 10~ 3 -10 3 s. A cyclotron line at ~ 30 keV has been suggested by various authors but 
never detected with significance. The stellar wind of the optical companion is abnormally slow. The orbital period, initially reported 
to be 9.6 days, disappeared and a new periodicity of 19.25 days emerged. 

Aims. The main objective of our RXTE monitoring of 4U 2206+54 is to study the X-ray orbital variability of the spectral and timing 
parameters. The new long and uninterrupted RXTE observations allow us to search for long (~1 hr) pulsations for the first time. 
Methods. We divided the ~7-day observation into five intervals and obtained time-averaged energy spectra and power spectral density 
for each observation interval. We also searched for pulsations using various algorithms. 

Results. We have discovered 5560-s pulsations in the light curve of 4U 2206+54. Initially detected in RXTE data, these pulsations 
are also present in INTEGRAL and EXOSAT observations. The average X-ray luminosity in the energy range 2-10 keV is 1.5 x 10 35 
erg s _1 with a ratio F max /F m j n ss 5. This ratio implies an eccentricity of ~0.4, somewhat higher than previously suggested. The power 
spectrum is dominated by red noise that can be fitted with a single power law whose index and strength decrease with X-ray flux. 
The source also shows a soft excess at low energies. If the soft excess is modelled with a blackbody component, then the size and 
temperature of the emitting region agrees with its interpretation in terms of a hot spot on the neutron star surface. 
Conclusions. The discovery of X-ray pulsations in 4U 2206+54 confirms the neutron star nature of the compact companion and 
definitively rules out the presence of a black hole. The source displays variability on time scales of days, presumably due to changes 
in the mass accretion rate as the neutron star moves around the optical companion in a moderately eccentric orbit. If current models 
for the spin evolution in X-ray pulsars are correct, then the magnetic field of 4U 2206+54 at birth must have been B > 10' 4 G. 

Key words, stars: individual: 4U 2204+54, BD+53°2790 - X-rays: binaries - stars: neutron - stars: binaries close -stars: emission 
line, Be 



1. Introduction 

The X-ray source 4U 2206+54 is one of the most enigmatic 
high-mass X-ray binaries (HMXB). Although it has been stud- 
ied by numerous ground- and space-based observatories, many 
questions about the nature of its variability patterns remain unan- 
swered. I n the optical and UV band s, the emitting spectrum is 
complex (Negueruela & Reig 200l|) an d far from that of a Be 
star. Initially classified as a BO-2 star dSteiner et al.l fl984), the 
optical counterpart defies any standard spectral analysis. The ab- 
sence of any obvious long-term trends in the evolution of the 
Ha line spectral parameters and the lack of any correlation be- 
tween the Ha line equivalent width and infrared magnitudes and 
colours rules out a Be classification. However, the strength and 
shape of the Ha line do not resemble those seen in supergiant 
stars either. Although the presence of an 09.5V st ar is virtu- 
ally s ecured, the He abundance is abnormally strong (Blay et al. 
2006). 
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Without a circumstellar disc around the donor, the material 
needed for accretion and production of X-rays must come from 
the stellar wind. However, the expected X-ray luminosity us- 
ing the typical wind velocities of an 09.5V star is about two 
orders of magnitude lower than the observed X-ray luminosity 
L x ~ 10 35 er g s~'. The distance to 4U 2206+54 is estimated 
to be 2.6 kpc dBlav et al.ll2006l) . iRibo et all d2006l) found that 
the wind terminal velocity of 4U 2206+54 is abnormally slow 
(~350 km s~'). With such a low velocity, an eccentric orbit and 
using the Bondi-Hoyle formalism it is possible to reproduce the 
average X-ray luminosity of the system as well as the orbital 
variability seen in the RXTE/ ASM light curve between 1996 and 
2005. 4U 2206+54 is the only permanent wind-fed HMXB with 
a main-sequence donor dRibo et al.l l2006). 

Equally enigmatic is the nature of the compact companion, 
which has not been settled beyond doubt. Although a neutron 
star is the most likely scenario, the low X-ray luminosity, com- 
pared to other persistent HMXB, has been used as an argument 
to suggest the presenc e of a white dwarf dSaraswat & Apparaol 
Il992i ICorbet & Peeldl200lh , while the lack of pulsations and 
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Fig. 1. Light curve of 4U 2206+54 covering the entire interval 
of the RXTE PCA observations. The data points represent the 
average of each observation interval as given in Table Q] Errors 
are less than the size of the points. Dashed lines separate the 
different intervals for orbital variability analysis (see text). The 
numbers on top of the figure are the orbital phase for Pi = 9.56 
d and Pj = 19. 12 d. Reference time for phase calculation is MJD 
53567.74, which corresponds to m aximum flux in the folded 
ASM light curve dCorbet et al.ll2007l) . 

Table 1. Log of the RXTE observations. 
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01-02-00 
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3776 


7.3 
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01-02-040 
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6.2 
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01-02-04 
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1984 
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a: Background subtracted 4-30 keV, PCU2 only 



the similarities with LS 5039 do not allow us to rule out a black 
hole companion dNegueruela & Reig 2001 ) - the most recent re- 
sults show, however, that LS 5039 is probably a non-accreting 
neutron star dRibo et al. 2008). 
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Fig. 2. Portion of ~ 6 hours of the light curves of 4U 2206+54 
corresponding to each one of the observation intervals. The axis 
scales were left the same in all panels to allow easy comparison 
of the variability amplitude. The bin size is 10 seconds. Intervals 
C and D are more affected by observational gaps. 



Broad-band X-ray observations favour the presence of 
a neutron star dTorreion et alj 120041) . When modelled with 
Comptonization models, the resulting seed photons have kT > 1 
keV, while the bolometric luminosity is low (Lx ~ 10 35 erg s _1 ). 
The only way to reconcile these two results is by invoking a 
small emi tting area (< 2 km ), like that of a polar cap in a neu- 
tron star (Masetti et al. 2004), also ruling out the presence of an 
accretion disc ( Torreion e t al. 2004). The presence of a magnetic 
field with B ~ 2 x 10 12 Gis implied by the possible detection 
of a cyclotron resonance scattering feature at -30 keV. Although 
the statistical significance of this feature is only marginal, it has 
been observed by three diff erent observatories, namely RXTE, 
BeppoSAX and INTEG RAL dTorreion etalj|2004t iMasetti et al.1 
l2^MlBlavlt^ni2005l) . 

Another debated questio n is the orbital period of the sys- 
tem. ICorbet & Peeld d200ll) reported an orbital period of 9.6 
days. This value was based on observations made with the RXTE 
ASM instrumen t over 5 years. Thi s orbital period was refined to 
9.5591 days bv lRibo etalj d2006) using 9 years wor th of data 
of the same instrument. However. ICorbet et al.l d2007l) could not 
find the 9.6-day modulation in the light curves of SWIFT BAT 
and RXTE ASM for observations carried out from 2004 to 2005. 
Instead, they found a strong modulation at 19.25 d, that is, con- 
sistent with twice the 9.6-day period (or 19.12 d). 

In this paper we investigate the X-ray variability of 
4U 2206+54 on time scales of hours to days. The main objec- 
tive is to search for ~l-hr pulsations and study possible spectral 
and timing orbital variability. 



2. Observations 

4U 2206+54 was observed by the Rossi X-ray Timing Explorer 
(RXTE) during the interval 15-22 May 2007 (JD 2454235. 1628- 
2454242.1991) for a total on-source time of 136.8 ks. 
Observation time with RXTE was awarded in response to the 
All call for proposals and corresponds to the RXTE proposal 
P92068. Light curves, spectra and response matrices were ex- 
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Fig. 3. Hardness-intensity diagram of 4U 2206+54 for various 
X-ray bands. Count rates correspond to PCU2. Bin size is 512 s. 



tracted using the tools available in the RXTE subpackage of ver- 
sion v.6.4 of the FTOOLS software package. 

RXTE is designed to facilitate the study of time variability in 
the emission of X-ray sources with moderate spectral resolution. 
The pay-load consists of a Proportional Counter Array (PCA), 
which is sensitive to X-rays in the 2-60 keV energy range and 
has a total effective area of 6500 cm 2 , shar ed by five PCU or 
proportional counter units ( Jahod a~et al.lll996h . PCA data can be 
collected and telemetered to the ground in many different ways 
depending on the intensity of the source and the spectral and 
timing resolution desired. In this work we used the two stan- 
dard modes: Standard! provides 0.125-s resolution and no en- 
ergy resolution; in the Standard2 configuration data are accu- 
mulated every 16 seconds in 129 channels. The High Energy 
X-ray Timing Experiment (HEXTE) comprises two indepen- 
dent clusters of detectors with a total collecting area of 1400 
cm 2 , covering the energy range 15-250 keV. The spectral res- 
olution is ~\5% at 60 keV energy range. During normal oper- 
ation, the tw o clusters rock on an d off target to collect back- 
ground data (Roths child et al.|[l998l) . Due to malfunction of the 
rocking mechanism, cluster A was left permanently on the on- 
source position for the entire duration of the observations. Thus, 
no background light curve could be obtained. Data presented in 
this work correspond to cluster B only. The third instrument on 
board RXTE is the All-Sky Monitor (ASM),which produce daily 
flux ave rages in the energ y range 1.3-12.1 keV at 30 mCrab sen- 
sitivity dLevine et alJll996l) . 

The log of the RXTE observations is shown in TableQ] Good 
time intervals were defined when the pointing of the satellite was 
stable (< 0.02° from the source), the elevation above 8° and far 
away from the South Atlantic Anomaly. Due to RXTE's, low- 
Earth orbit, the data consist of a number of contiguous data in- 
tervals (typically 1 hr long) interspersed with observational gaps 
produced by Earth occultations of the source and passages of the 
satellite through the South Atlantic Anomaly. In general, the vis- 
ibility of a target depends on its position with respect to the plane 
of RXTE's, orbit around the Earth. Objects that are within a crit- 
ical angle from the normal to the orbit will be visible for longer 
periods of time. In fact, if the object is close enough to the nor- 
mal then it would lie in the so-called Continuous Viewing Zone 
(CVZ). Since RXTE orbits relatively close to the Earth's equa- 



tor, the regions of high visibility and the CVZ are close to the 
north and south pole of the Earth. 4U 2206+54 is at high ecliptic 
latitude (about 58 degrees), hence it has a relatively high vis- 
ibility. The observations reported here are considerably longer 
than commonly reported for other sources. There are two long 
(~ 6 hr) uninterrupted observations. These longer observations, 
together with the monitoring of 4U 2206+54 over 7 days, allow 
us to give new insights into the properties of the X-ray emission 
in 4U 2206+54. In particular, they allow us to search for peri- 
odicities on timescales of ~ 1 hr. To investigate potential orbital 
variability, the RXTE observations were divided into five inter- 
vals with roughly the same on-source time, ~20-30 ks (Fig.Q]). 

In addition to the RXTE pointing observations, we also anal- 
ysed archived data from the INTEGRAL and EXOSAT mis- 
sions. The International Gamma-ray Astrophysics Laboratory 
(INTEGRAL) observations were made between 16-19 December 
2006 during revolution 510. INTEGRAL carries two main 
gamma-ray instruments: the high angular resolution Imager 
(IBIS) and the high-energy resolution spectrometer (SPI). In our 
study we only used data from the upper detector system of IBIS, 
called the Soft Gamma-ray Imager (ISGRI). ISGRI is sensitive 
to photons in the energy range 20 keV to 1 MeV, has an angu- 
lar resolution of 12' and a source location accuracy that depends 
on the significance of the detection: for a source at the limit of 
detection the source location accuracy is abo ut 4-5' while for a 
strong source it can go down to about 0.5' ( Gros et aL| l2003). 
The energy resolution is 7% at 100 keV dUberfini et all 120031) . 
4U 2206+54 was detected in a total of 19 INTEGRAL science 
windows of around ~2.3 ks exposure each. In all these obser- 
vations the satellite followed a 5 x 5 dithering pattern. We only 
analysed images in which the source fell in the fully coded field 
of view (9° x 9°). Data reduction was done with the Offline 
Scien tific Analysis (OSA) software version 7 dCourvoisier et al.l 
120031 

The European Space Agency's X-ray Observatory, EXOSAT, 
was operational from May 1983 to April 1986. The pay-load 
consisted of 2 Wolter Type I grazing incidence Low Energy 
(LE; 0.05-2 keV) Imaging Telescopes, a Medium Energy (ME) 
Proportional Counter and a Gas Scintillation (GS) Proportional 
Counter. The ME detectors covered the energy range 1-50 
keV and had a field of view of 45 arcmin and an effective 
area 1600 cm- dTurneretalJU981l) . For our study we retrieved 
the ME detector standard products available in the NASA's 
HEASARC archive. EXOSAT observed 4U 2206+54 with the 
Medium Energy (ME) proportional counter in three occasions: 8 
August 1983, 7 December 1984 and 27 June 1985 for a total of 
9.8 ks, 9.3 ks and 6.3 ks, respectively. 

3. Timing analysis 

We have performed a timing analysis on two time scales: a) min- 
utes to hours to search for slow pulsations and b) days to inves- 
tigate orbital variability. Figure Q] shows the entire RXTE/PCA 
light curve. Each point represents the average of one observation 
interval (see TableQ]). On short time scales the X-ray light curve 
of 4U 2206+54 is dominated by erratic flaring, with changes in 
the X-ray intensity by a factor 3-5 over a few minutes (Fig. [2j. 
Despite the different wind paramete rs and optical companion be- 
lieved to be present in 4U 2206+54 dBlav et al.ll2006tlRib6 et al.l 
2006), the light curves are very similar to those of X-ray sources 
in which a neutron star is fed by direct accretion from the wind 
of a supergiant, such as Vela X-l. This supports the idea that ac- 
cretion in this source is not mediated by an accretion disc, but 
proceeds directly from the wind. Figure [2] shows a portion of the 
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light curve for each interval. The Y-axis spans the same intensity 
range in all panels to show the large variability of 4U 2206+54 
on short timescales. 

Figure|3]shows a colour-intensity diagram where three hard- 
ness ratios are plotted as a function of 4-30 keV PCU2 count 
rate. Each point represents an average over 512 seconds. The 
source becomes harder as the count rate increases. A fit to a con- 
stant function does not give good fits to the top two panels, with 
reduced^ of 8.5 and 5.1, indicating that the variation is statis- 
tically significant. In contrast, the data of the bottom panel give 
X; — 1 .4. The correlation coefficients to a linear regression fits 
in the three panels are, from top to bottom, 0.6, 0.7 and -0.2. 

3.1. Search for long-term pulsations 

In spite of intensive searches, X-ray pulsations have not been 
detected in 4U 2206+54. Searches for X-ray coherent emis- 
sion have been performed with EXOSAT dCorbet & Peelell200lL 
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but se e ISaraswat & App arao 1992 |), RXTE (Negueru ela & Reigl 
l200lh iTorreion et all 12004 ICorbet et all 120071). BevvoSAX 
iMasetti et alJ l2004t rTorreion et all 120041) and INTEGRAL 



Blay et al. 2005). These studies demonstrate the lack of X-ray 
pulsations on timescales from ~ 1 ms to ~ 1 hour. The possible 
detect ion of 392-s pulsations rep orted by Saraswat & Apparaol 
(i!992h has never been confirmed. Corbet & Peele (2001) reanal- 
ysed t he EXOSAT observations used by ISaraswat & Ap parao 
(1992) and did not find any evidence of any pulsat ions longer 
than 2 s. They argue that Saraswat & Apparaol (1 1992b used a pro- 
cedure not well suited to search for coherent periodic pulses. 

The only timescales that have not been investigated in detail 
due to either technical constraints (low-Earth orbits) or lack of 
long uninterrupted observations is from ~ 1 hr to ~ 1 day. The 
long coverage of our RXTE observations allows us to investigate 
the higher frequency part of this unexplored interval. 

We searched for periodic signals in the light curves us- 
ing various standard techniques. Initially we applied an FFT to 
the two longest observation intervals (see Table [TJ. However, 
in order to be able to use the entire data set, techniques suit- 
able for handling missing data (gaps) i n the light curves ar e 
needed. Thus we employed PERIOD04 dLenz & Breger 2005 
the Phase Dispersion Minim isation (PDM. IStellingwerf|[l97.. 
and the CLEAN algorithm dRoberts et alj|1987l) . A periodic sig- 
nal around ~5560 s is clearly detected with all methods. The 
PDM and CLEAN periodograms obtained between 100 and 
10000 s for the RXTE/PCA light curve binned at 16 s intervals 
are shown in the top panel of Fig. [4] In the case of PDM there is 
a strong beating with ~400 s separation, and two equally signif- 
icant minima occur at 5555 ± 10 s and 5608 ± 10 s. In the case 
of CLEAN there is a clear maximum only at 5559 ±3 s, 

To further study the periodic signal detected with the whole 
data set, we have analysed the RXTE/PCA light curves of inter- 
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Fig. 6. Power spectra of 4U 2206+54 at different instances of the 
observations. See Table|2]for a definition of the intervals. 



vals A+B and intervals C+D+E separately, corresponding to a 
high and low flux states of the source on an orbital timescale 
(Fig.Q}. The results show that the periodic signal is detected in 
both data sets, clearly indicating that it is independent of the X- 
ray flux of the source. In particular, for A+B PDM and CLEAN 
provide signals at 5550 + 50 s and 5562 + 3 s, respectively. For 
C+D+E we obtain 5560 + 30 s and 5559 ± 3 s, respectively. 
The probability that the detected period is due to random fluctu- 
ations is less than 1 0~ 8 , according to th e statistical properties of 
the PDM algorithm dStellingwerf 19781) . 

We note that a strong signal is also detected around 
815 000+15000 s (9.4 + 0.2 d) with the CLEAN algorithm in 
the RXTEfPCA light curve. Although the available database is 
limited to just 7.0 d, it is interesting to note that a clear hint of 
the possible 9.6 d (or of a part of the 19.25 d) orbital period 
appears in the RXTEfPCA data set of 4U 2206+54. The proba- 
bility that there is not such periodicity or that the real period is 
different to 815 000 s is less than 0.01. 

The value of 5560 s is very close to RXTE's orbital period - 
RXTE follows a low-Earth circular orbit at an altitude of 580 
km, corresponding to an orbital period of about 90 minutes. 
Note that the RXTE's orbital period is continuously changing 
due to perturbations by the Earth and Moon and atmospheric 
drag. At the time of the observations the orbital period laid in 
the range 5667.3-5669.1 s. In order to rule out an instrumen- 
tal origin of the modulation, we analysed archived data from the 
INTEGRAL/ISGR1 and EXOSAT/ME missions and instruments. 
The orbit of these two observatories is highly eccentric with a 
period of 72 hours and 90 hours, respectively. The time bin was 
500 s for the ISGRI light curve and 30 s for the ME light curve. 

The timing analysis on the INTEGRAL (20^0 keV) and 
EXOSAT (0.8-9 keV) light curves gives periodic signals at 
5570 ± 20 s and 5525 ± 30 s, respectively (Fig.[4]i, which clearly 
confirms the presence of the 1 .54-hr periodic signal found with 
the RXTE/PCA data. 

Figure[5]shows the pulse profile corresponding to each satel- 
lite. The pulsed fraction is -50% for the RXTE and INTEGRAL 
light curves and 35% for the EXOSAT data. 

HEXTE data were not used in the periodicity analysis be- 
cause of the problems with cluster A (see Sect. 2) and the low 
S/N of the cluster B light curves. Also, note that the energy range 
covered by HEXTE is partially covered by ISGRl/INTEGRAL. 
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3.2. Orbital variability 

To investigate the longer-term variability, the observations were 
divided into five intervals with roughly the same on-source time. 
The average intensity decreased throughout the time covered by 
the observations as can be seen in Fig. [TJ The numbers on top 
of this figure are the orbital phas e assuming Pi = 9.56 d and 
twice that period P 2 = 19.12 d (ICorbet et al.ll2007l) . The refer- 
ence time is the epoch of maximu m flux as given by the ASM, 
MJD 53567.74 (ICorbet et al.ll2007l) . The PCA light curve shows 
maximum flux at phase ~0, in agreement with t he modulation 
of the folded ASM light curve (ICorbet et al.ll2007l) . There is an- 
other less significant peak at about phase 0.5, if P or b = 9.56 d. If 
the orbital period were P 01 -b = 19.12 d, then we would expect t o 
see the maximum flux around phase 0.125 dCorbet et al.ll2007l) . 
which is not seen. 

In order to investigate the evolution of the power spectral pa- 
rameters of the noise components as the source moves around 
the orbit, we produced one power spectrum for each interval. 
The power spectra were obtained by dividing the 2-20 keV PCA 
2~ 7 -s binned light curve of the entire observation into 1024-s 
segments and an FFT obtained for each segment. The contribu- 
tion by the photon counting noise was computed and subtracted 
from each power spectrum. The power spectra were normalized 
such th at their integral gives the squared fractional rms vari- 
ability dBelloni & Hasingerlll990t iMivamoto etai1ll991l) . The 
power spectra of the same interval were averaged to produce the 
final power spectrum. A logarithmic rebinning was also applied 
to reduce the noise at high frequencies. Figure|6]shows the power 
spectra for the five intervals. The results of the fitting procedure 
are shown in Table [2] 

A single power-law component provides good fits for all in- 
tervals. The power-law index and the fractional rms amplitude 
change with X-ray intensity. The high-intensity intervals (A and 
B) show a steeper power spectrum and a larger amplitude of vari- 
ability. The low-intensity intervals (D and E) display a natter 
power law and are less variable. Interval C shows intermediate 
values between these two cases. 
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Table 2. Results of the power-law fits to the power spectra. 





MJD 




Orbital 


Number 


r 


rms (%) 


* r 2 /d0f 






c/s 


phase* 


of PSD 




(0.001-1 Hz) 




A 


54235.16-54236.11 


29 


0.87/0.93 


21 


1.87±0.04 


37 


1.3/44 


B 


54236.86-54238.00 


27 


0.06/0.03 


26 


1.77+0.04 


34 


1.1/44 


C 


54238.00-54239.06 


12 


0.17/0.09 


19 


1.55+0.05 


30 


1.9/43 


D 


54239.09-54241.03 


6 


0.33/0.17 


24 


i c+0.09 
1 - J -0.13 


28 


1.3/44 


E 


54241.05-54242.20 


9 


0.50/0.25 


15 


1.5+0.1 


23 


1.2/44 



a: Background subtracted 2-60 keV, PCU2 only 

b: At mid-point for P orb = 9.56 d and P orb = 19.12 d, respectively 



4. Spectral analysis 

In this section we present the results of the energy spectral anal- 
ysis. We first present the PCA+HEXTE time-average spectrum 
for the entire observation. Then, we search for orbital variability 
in the spectral parameters. 



4.1. Time-averaged spectrum 

The 3-30 keV spectrum of 4U 2206+54 is well described by 
an absorbed, Ah = 1.3 + 0.5 x 10 22 cirT 2 , power law with 
photon index 1.3+0.1 modified by a cutoff at 6.3+0.4 keV and 
folding energy 14+2 keV , in agreement with t hose found by 
iNegueruela & R eig (2001). However, as shown in Torreio rTet alj 
d2004l) . this model does not describe properly the spectrum when 
data beyond 30 keV are taken into account (y 2 =1.44 for 53 
d.o.f.). The joined time-average PCA+HEXTE spectrum (Fig.[7]l 
is, overall, well described by thermal Comptonization (comptt 
in xspec notation). This model requires further a blackbody emis- 
sion (bb) component to account for an excess at low energies. 
The same is true when dynamical Comptonization (bmc) re- 
places thermal Comptonization. The parameters for these two 
composite models are presented in Table [3] kT\,b is the temper- 
ature of the blackbody component that accounts for the soft ex- 
cess; kTo and kT co \ are the temperatures of the seed photons, 
which are subsequently Comptonized by the high temperature 
plasma, in the thermal and dynamical Comptonization models 
respectively; kT e and t are the temperature and optical depth 
of the high temperature plasma, in the thermal Comptonization 
model; a is an spectral index related to the efficiency of the 
Comptonization process (the lower the value of a, the more ef- 
ficient is the Comptonization); / is the illumination factor, i.e., 
the ratio of the number of multiply scattered photons to the total 
number of seed photons (f » 1 implies that the spectrum is fully 
Comptonized). No iron line is required. If the absorption column 
is let free to vary, then the best-fit value is higher than in previ- 
ous observations. Fixing it to a lower value increases slightly the 
temperature of the blackbody but it is always around 0.5 keV 
(for example fixing Ah to 0.3 x 10 22 cirT 2 , the pure interstellar 
value, gives kT^ = 0.56 and sets the^ 2 = 1). Since RXTE can- 
not constrain the column density due to the lack of sensitiveness 
below 3 keV we fixed Ah to the value repo rted from BeppoSAX 
observations, name ly A H = 0.9xl0 22 cm ~ 2 dMasetti et al.l2 004). 
Applying Eq. 1 in Torr eion et alj d2004l) to the data in Table [3] 
we determine the radius of the soft photon source which would 
produce the observed luminosity^ to be Rw = 1.5 km. This is 
only consistent with a hot spot on the neutron star surface. 



Table 3. Parameters for the Comptonization models. Error bars 
are 90% confidence level. 



1 This formula gives the proper factor for a circular spot of area nR^. 
Of course, a sphere would require a radius R sph = R spot /2 to produce the 
same L x assuming isotropic emission. 
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(a) 

{b) 3-60 keV unabsorbed flux in units of 10~ 10 erg s~' cm 



Hydrogen column density in units of 10 22 cm 2 (fixed) 




Pulse phase 

Fig. 8. X-ray flux and photon index as a function of the spin 
phase. 



The possi ble cyclotron line around 30 keV reported by vari- 
ous a uthors dTorreion et al.ll2004trMasetti et alJl2004tlBrav et al] 
120051) is not seen in the present data set. We did not find sig- 
nificant evidence for the iron emission line at ~6.4 keV either. 
We estimate an upper limit on the equivalent width of a narrow 
line (FWHM=0.1 keV) around 6.4 keV of 25 eV at 90% confi- 
dence level. The average X-ray luminosity in the energy range 
3-30 keV is 2 x 10 35 erg s~', assuming a distance of 2.6 kpc 
dBlav et al.ll2006l) . 
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Table 4. Results of the exponentially cutoff power law fits to the energy spectra. In the bottom table A^h was kept fixed. 
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(a) : At mid-point for P OI +, = 9.56 d and 2.P OI +, = 19.12 d, respectively 

(b) : Hydrogen column density in units of 10 22 cirT 2 

(c) : Unabsorbed flux in units of 10~ 10 erg s -1 cirr 2 



4.2. Pulse-phase spectroscopy 

In order to investigate how the X-ray spectrum changes as a 
function of the spin of the magnetized neutron star, the pulse 
profile corresponding to intervals A and B was divided into eight 
equally spaced bins and a spectrum in the energy range 3-20 keV 
was obtained for each bin. After subtracting the underlying per- 
sistent emission, each spectrum was satisfactorily fitted with an 
absorbed power law. The column density was fixed to 1 .3 x 10 22 
cm -2 , i.e., the value found in the time-averaged spectrum (Sect. 
4.1), although the use of A^h = 0.9 x 10 22 cirr 2 did not change 
the results significantly. Figure[8]shows the variation of the 3-20 
keV X-ray flux and the power-law photon index as a function 
of the pulse phase. Phase zero corresponds to MJD 54235.169. 
Marginal evidence for variability in the power-law photon in- 
dex is found. The spectrum becomes softer as the pulse phase 
increases, with the photon index reaching a maximum around 
phase 0.9. The hardness of the spectrum does not seem to de- 
pend on the pulse flux, but the photon index appears to lag the 
pulse flux by about 1200 s. However, the uncertainty in the pho- 
ton index values is too large to draw a firm conclusion. 

4.3. Orbital spectral variability 

To investigate possible orbital spectral variability, we obtained a 
time-averaged spectrum for each observation interval. The 3-30 
keV band spectrum was fitted to an absorbed power law modified 
at high energies by an exponential cutoff. Given the high low- 
energy limit, the hydrogen column A^h is not well constrained. 
Although the observations that correspond to a lower flux have 
lower column densities, the value of A^h is consistent with no 
variability, within the errors. Thus, we also tried fits with a fixed 
A^h = 0.9 x 10 22 cirr 2 (lower part of Table |4j. In either case, 
a variation in the power-law index is apparent. High-flux states 
show harder spectra. 

5. Discussion 

We have reported the analysis of the monitoring campaign of 
the high-mass X-ray binary 4U 2206+54 that took place in May 
2007 with RXTE. The observations span over ~ 7 days, which 
represent a large fraction of the 9.6-day period. The long and 
uninterrupted observations allow us for the first time to i) ex- 
tend the search for X-ray pulsations to time scales of about one 
hour, ii) investigate the X-ray orbital variability and Hi) give 



new insights into its energy emission properties. In support of 
the results of these observations we also analysed data from 
INTEGRAL/ISGR1 and EXOSAT/ME. 

5.1. X-ray pulsations 

The proximity of 4U 2206+54 to RXTE's Continuous Viewing 
Zone allowed the acquisition of long and uninterrupted observa- 
tions, and the discovery of slow X-ray pulsations (P sp i n = 5559 + 
5 s). With the discovery of X-ray pulsations we have solved one 
of the questions that remained unanswered, namely, the nature 
of the compact companion. The vast majority of HMXBs har- 
bour X-ray pulsars. These are believed to be young neutron stars 
with relatively strong magnetic fields (B ~ 10 12 G). Among the 
handful of HMXBs not displaying X-ray pulsations, three show 
the typical characteristics of persistent black-hole systems. Of 
these, two are in the LMC and only one (Cyg X-l) in the Galaxy. 
The nature of the other three HMXBs in which pulsations have 
not been discovered in spite of intensive searches is under de- 
bate, although two of them may also cont ain black hole com- 
panions: LS I +61 303/2E 236.6+6101 (ICasares et al.ll2005al 
but see iDhawan et all |2006|) and LS 5039/RX J1826.2-1450 
(ICasares et al.ll2005bl but see iRibo et all [20081) . If confirmed it 
would explain the lack of pulsations in these systems. Geometric 
effects (very low inclination or a very small angle between the 
spin and rotation axises) have been invoked to explain the lack of 
pulsation s in 4U 1700-37, which might conta in a massive neu- 
tron star (IClark et al.ll2002t lAbubekerovll2004 . The lack of pul- 
sations was more difficult to explain in 4U 2206+54, as previous 
observations indicated that this source likely hosts a canonical 
neutron st ar accreting from the wind of a main-sequence early- 
type star dTorreion et ail 12004 iMasetti et all 120041; iBlav et alj 
120051 120061: IRibo et al.ll2006l) . 

With the period detected here, 4U 2206+54 is the third (if 
the nature of IGR J16358-4726 as a HMXB is finally con- 
firmed) slowest HMXB pulsar, after 2S 0114+650 and IGR 
J16358^4726, and the first with a main-sequence c ompanion. 
2S 0114+650 is a 2.78-h pulsar jFinley et al.lll994l) ac creting 
matter from a supergiant Bl companion dReig et al Il996l). while 
IGR J 16358- 4726 has been clas sified as a sgBlel (IChaty et al.l 
120081 but seelNespoli et al.1 12008) with a 1.6-hr rotating neutron 
star dPatel et al . 2004)0 



The slowest X-ray pulsar is 4U 1 954+319 with P. 



5 hr 



(Matt ana et"aIH 2006: Corbet et al. 2008). However, this is a symbiotic 
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Numerical calculations suggest that there is no signifi- 
cant angular momentum transfe r onto the ne utron star from 
the wind of supergiant systems (Ruff era U999h . Even a lower- 
velocity wind like the one believed to be present in 4U 2206+54 
(iRiboetalJ 120061) seems to be insufficient to exert a spin-up 
torque onto the neutron star. Since the transfer of angular mo- 
mentum during the accretion phase is so inefficient, a spin period 
of ~ 5500 s must have been attained in a previous evolutionary 
phas e. The neutron star's spin evolution in a close binary sys- 
tem dDavies et alJl979HDavies & Pringldl98Tl:lDai et al.ll2 006) 
is divided into three phases, of which the accretion phase is the 
more recent one. The neutron star is born as a millisecond pul- 
sar (pulsar phase) and spin down ocurs by magnetic dipole ra- 
diation until P S pi n ~ 1 s. Then the neutron star enters the pro- 
peller phase, in which the neutron star spins down because ac- 
cretion is centrifugally inhibited. This phase continues until the 
spin period reaches the so- called equilibrium period give n by 
dLi & van den Heuvelll 19991 but see lDavies & Pringldl 19811) 

P eq » 20< 7 M,f X 8/7 M-.4 /7 s 

where B = 10 12 Pi2 G is the neutron star's dipolar magnetic 
field strength, M = 10 15 Mi5 g s is the mass accretion rate, 
R = 10 6 /?6 cm is the radius of the neutron star, and M = 
1 AM i .4 M Q is the mass of the neutron star . Adopting the c anon- 
ical values R 6 = Mi 4 = 1 and B 12 = 3.6 dBlav et alj|2005l) and 
Mis ~ L X R/GM = 0.5 for L x = 10 35 erg s"\ the 4U 2206+54 
equilibrium period results P eq ~ 80 s, far lower than the ob- 
served period. We conclude then that the current spin period of 
4U 2206+54 cannot be explained in evolution ary terms. A way 
out is provided bv lLi & van den Heuvell dl999h . who explain the 
long spin period of 2S 01 14+65 assuming that the neutron star 
was born as a magnetar, that is, with B > 10 14 G. Indeed, the pro- 
peller effect can spin down the neutron star to 5 x 10 3 s keeping 
the same values as above but taking B\2 two orders of magnitude 
higher. 

5.2. Temporal variability 

The range in X-ray luminosity (3-30 keV) implied from the ob- 
servations is 0.6-3 .1 xlO 35 erg s" 1 , assuming a distance to the 
source of 2.6 kpc dBlav et al.ll2006l) . That is, 4U 2206+54 dis- 
plays an amplitude of variability in X-ray flux of ab out 5. This 
amplit ude is larger than the factor of 2 reported by Rib o et al] 
d2006l) from ASM data, i.e., for X-rays i n the 1.3-12 keV band 
but consistent with the observations by ICorbet & Peelel ([2001, 
see their Fig. 4). Note that the results from the ASM are highly 
dependent on the minimum flux, which for the ASM is uncer- 
tain. In the simplified Bondi-Hoyle accretion model a maximum 
to minimum orbital X-ray flux ratio of 5 can be explained as- 
suming that the orbital variability is due to changes in the mass 
accretion rate. In this scenario, an eccen tricity of ~ 0.4 is n eeded, 
somehow higher than that suggested by Rib o et al.l (120061) . 

In addition to the X-ray flux, the power and energy spectral 
parameters also display variability on time scales of days. If the 
orbital period of 9.6 days is assumed to be correct, then it seems 
natural to attribute this variability to the motion of the source 
through the orbit. Note however, that our results are based on 
less than one orbit. The power spectra of 4U 2206+54 are well 
described by a single power law, whose index and rms ampli- 
tude decrease as the count rate decreases (Table [2]). There is no 

X-ray binary with a neutron star c omponent accreting from a low-mass 
M-type giant jMasetti et alj2006bh . 



sign of a break at low frequencies (flat-topped noise) in the fre- 
quency range covered by the power spectra. Red noise dominates 
the power spectrum down to frequencies as low as 1 mHz. The 
decrease in the amplitude of the X-ray variability (rms) with de- 
creasing X-ray flux seems to indicate a large scale change in the 
structure of the wind. One possible explanation is that the wind 
was highly inhomogeneous during intervals A, B and C, during 
periastron passage, as usually observed in systems with super- 
giant companions, and became smoother further out, in intervals 
C and D. Assuming a mass of 18M for the 09.5 V primary and 
the canonical neutron star mass of 1.4M Q for the accretor, the 
change in the wind structure must have occured in between in- 
tervals C and D, which in terms of radial distance corresponds to 
-5.7/?* for Pi = 9.6 d and -6.95 P, for P 2 = 19.25 d. The detec- 
tion of X-ray emission in intervals D and E in 4U 2206+54 con- 
trasts with what is observed in classical supergiant systems for 
which the X-ray emission is hardly detected when the compact 
objec t is located beyond 3 P* dNegueruela et al.l2 008; Blav et all 
120081 

The 3-30 keV spectrum of 4U 2206+54 can be described 
with a power law with T = 1.5 and a high-energy cutoff with 
Pcut = 6.7 keV (see Sect. 4.1). The power-law index increases 
as the X-ray flux decreases (Table 0). The high-flux intervals 
(A and B) show a harder spectrum than the low-flux states (C, 
D and E). This behaviour is the opposite to that found in other 
X-ray sources. The explanation could be, again, in the low X- 
ray brightness of this source. Within the framework of bulk mo- 
tion Comptonization, that is, assuming that Comptonization oc- 
curs in the accretion flow near the neutron star surface (accretion 
column), the surface of the neutron star acts as source of feed- 
back against the infalling material: when the accretion increases, 
so does the X-ray flux. The ram pressure exerted by the X-ray 
source tends to hamper accretion and therefore complicates the 
Comptonization process. In the very bright X-ray sources, the 
feedback decreases the efficiency of Comptonization and pro- 
duces a high-soft low-hard behaviour. In 4U 2206+54 however, 
with its main sequence donor, the X-ray source is never pushed 
to the limit where the X-rays can hamper the accretion process. 
With no feedback from the surface of the neutron star, the source 
becomes harder when accreting material becomes available, pre- 
sumably at the periastron passage. 

One of the controversial results of 4U 2206+54 in the 
X-ray band is the shape of the hardness-intensity diagram. 
Negueruela& Reig (2001) found a correlation between the X- 
ray colours and the X-ray intensity, in that the X-ray spectrum 
becomes harder as the flux increases. Figure [3] confirms these 
results. This figure shows the hardness-intensity diagram for 
three different ratios. In all three cases the change in intensity 
is approxima t ely the same, namely, a factor of 10. In contrast, 
Mas etti et al.l d2004l) did not find any variability of the hard- 
ness ratio with intensity i n a B eppoSAX observation on 1998 
November. iMasetti et al.l (120041) suggested that the hardness- 
intensity dependence might become evid ent at higher lumi- 
nositie s, as the X-ray luminosity during the iNegueruela & Reigl 
(2001) obse rvations was about one order of magnitude higher 
than that of Masetti ~et al.l d2004l) . To investigate this possibility 
we obtained hardness-intensity diagrams for interval B and D, 
that is, when the X-ray flux was maximum and minimum, re- 
spectively. We define the hardness ratios as HR1 =7-10 keV/4- 
7 keV and HP2 =10-15 keV/4-7 keV. Fitting HRl and HR2 as 
a function of 4-30 keV intensity to a constant does not give ac- 
ceptable fits for interval B. The fits to interval D are acceptable 
(in the sensed < 2). However, the reason for this acceptable fit 
is because the error bars are larger. As an example, a linear fit 
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Fig. 9. Variation of the characteristic temperature (upper panel) 
and radius (lower panel) with X-ray luminosity of the emitting 
region producing the soft excess, if modelled with a blackbody. 
Squares represent high-luminosity disc-fed supergiants (Cen X- 
3, SMC X-l and LMC X-4); filled circles are low-luminosity 
persistent Be/X-ray binaries (X-Per, A 0535+262 and LS I +61 
235); triangles correspond to the wind-fed supergiants (Vela X-l 
and 2S 0114+65); the open circle is SAX J2103.5+4545; stars 
mark the position of 4U 2206+54 for three different observations 
(see text for details). 



(y — ax + b) to the //7?2-intensity diagram of interval D reduces 
the x 2 from 419 (259 dof) to 290 (258 dot). An F-distribution test 
gives a chance improvement of only 0.18%. That is, the addition 
of a linear component is significant above a 3-cr level. In fact, ir- 
respective of the X-ray flux, as the intensity increases by a factor 
10, the hardness ratios HRl and HR2 increase by ~25%. In sum- 
mary, the hardness-intensity depe ndence becomes mo re evident 
at higher source luminosities (as Mas etti et al.l {2004) had sug- 
gested). However, this effect is most likely due to low statistics 
rather than the disappearance of the correlation. 



5.3. Soft component 

When a broader energy range is considered, the high-energy 
cutoff power-law model does not give acceptable fits and a 
soft excess at energies below 2 keV shows up in the residuals. 
More physical models involve Comptonization. Either thermal 
or bulk-motion Comptonization produce good fits to the X-ray 
energy continuum in the 3-80 keV band. In either case a black- 
body component with kTbb = 0.5 keV is required to account for 
the low-energy photons. 

The pre sence of a blackbody component supports the sug- 
gestion by iHickox et al.l (120041) that soft excess is common 
(probably ubiquitous) in X-ray pulsars. For high luminosity 
(Lx ~ 10 38 erg s _1 ) pulsars, the soft component stems from 
reprocessing of hard X-rays from the neutron star in the inner 
edge of an accretion disc. In less luminous systems, other mech- 
anisms, such as emission from diffuse gas or from the neutron 
star surface, can be at work. To get a rough estimate of the size of 
the soft radiation emission site, we compute the radius of the cir- 
cular blackbody spot that would produce the observed soft lumi- 
nosity, according to data in Table|3] We get a radius of Rbb ~ 1.5 
km. This value is too small to come from the entire surface of the 
neutron star, a hot corona, or an accretion disc. The lack of a disc 



in 4U 2206+54 was discussed in lTorrejon et al.l (|2004). This size 
is of the same order as that deduced previously from the temper- 
ature of the seed soft photons Ry/ a 1 .5 km. Then we tried to tie 
the temperature of both components during the fits. The result- 
ing statistics were significantly worse. It seems, therefore, that 
the spectrum comes from two sites, with the same characteris- 
tic size but different temperature. They could reflect a gradient 
in temperature from the hot spot (of radius /?w in the base of 
the accretion column) and the surface surrounding the spot (just 
outside the column) and into the relatively cold surface. 

When the soft excess is modelled with a blackbody compo- 
nent, the blackbody spectral parameters (temperature and emis- 
sion radius) strongly depend on the source luminosity, which in 
turn, depends on the accretion mechanism. In Fig. [9] we plot 
the temperature and radius of the emitting region allegedly caus- 
ing the soft excess. High-luminosity systems display lower tem- 
peratures and larger radii (open squares). In high-luminosity 
pulsars, such as SMC X-l, Cen X-3 and LMC X-4, accre- 
tion occurs via Roche lobe overflow and an accretion disc is 
formed. In these systems kT - 0.1 keV and j ?bb=100 km, typ- 
ically (see e.g. Table 5 in IHickox et al.ll2004l) . and the soft ex- 
cess is believed to be due to reprocessing of hard X-rays by 
optically thick accreting material. In cont rast, low-luminosit y 
systems (fill ed circles), such as X P er (ICoburn et al.l 1200 lb . 
A 0535+26 dMukheriee & Paull 12 005) and RX J0146.9+6121 
dLa Palombara & Mereghettil l2006h ~have kT = 1.2 keV and 
Rbb = 0.1 km. In these systems the neutron star captures ma- 
terial from the high-density wind from the circumstellar disc of 
the Be companion. The origin of the soft excess would be (part 
of) the neutron star surface. Note also that the high-luminosity 
pulsars mentioned above contain supergiant companions, while 
the low-luminosity systems harbour main-sequence stars. In be- 
tween these extreme cases one finds the wind-fed supergiant sys- 
tems (filled triangl es), such as Vela X -l, with kTbb = 0.2 keV 
and R bb = 70 km (IHickox et alj|2004l) and 2S 01 14+650, w ith 
kT bb = 0.3 keV and R bb = 200 km dMasetti et alj|2006al) . In 
these systems the thermal component responsible for the soft ex- 
cess is likely to be a clo ud of diffuse plasma around the neutron 
star (IHickox et al.ll2004l) . 

The star-like symbols in Fig. |9] represent thr e e different ob- 
serva tions of 4U 2206+54 dMasetti et alj 120041; iTorrejon et ail 
120041 and present work). 4U 2206+54 appears close to the Be/X- 
ray binaries, as expected if the soft excess arises from a hot spot 
on the neutron star surface. It is not clear how such a hot spot can 
be formed in 4U 2206+54 as accretion is also driven through a 
stellar wind dRibo e t al. 2006). Perhaps the much slower (and 
presumable denser) wind in 4U 2206+54 might make the accre- 
tion configuration be more similar to that of persistent Be/X-ray 
binaries. 

The open circle c orresponds to the peculiar Be/X-ray binary 
SAX J2103.5+4545 dlnam et alj|2004l) . This source is similar to 
4U 2206+54 in that the spin and orbital periods are typical of 
supergiant X-ray binaries, but the primary component of the bi- 
nary is a main-sequence star. In both systems the orbits have 
e = 0.4. The X-ray variability is, however, rather different. SAX 
J2 103. 5 +4545 presents X-ray bright states that last for a few 
months and extended (few years) X-ray faint states. During the 
bright states an acc retion disc is formed around the neutron star 
(Baykaletal. 2002) and a decretion circumstellar disc is formed 
around the Be star dReig et alj|2004l) . These discs are probably 
short lived and appea r during the high X-ray emission states 
onlv. lReig et"aIld2005l) found that SAX J2103.5+4545 was emit- 
ting X-rays even after the complete loss of the circumstellar disc. 
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Fig. 10. Pspm-Porb diagram. The two joined star-like symbols 
represent 4U 2206+54 for the two suggested orbital periods. The 
open circle corresponds to SAX J2103. 5+4545. 



hence further monitoring campaigns covering a larger fraction 
of the orbit are needed to pin down the correct orbital period. 
However, the fact that the maximum of the PCA light curve oc- 
curs at phase zero, i.e., in coincidence with the ASM maximum 
when folded onto the 9.6-day period and the detection of a strong 
signal at ~8 15000 s in the RXTE periodogram would favour the 
shorter period. The detection of a soft component in the energy 
spectrum has been interpreted as emission from the polar caps. 
The plot of the blackbody emitting radius and temperature as 
a function of X-ray luminosity constitutes a useful tool to dis- 
tinguish the type of accretion mechanism at work in the X-ray 
binaries that show a soft excess. 
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In this state. iReig et al.l (12005) argue that the X-rays are the result 
of wind-fed accretion, that is, as in 4U 2206+54. 

The position of these two sources in the P sp m-P or b diagram 
(Fig. [Tol l reinforces their similarity. They fall in the wind-fed 
supergiant region. Since th is diagram re flects the type of mass- 
loss/accretion mechanism ( Corbet] [l986l) . it is not surprising to 
find 4U 2206+54 and SAX J2103.5+4545 in this region. What 
makes them peculiar is that they both contain a main-sequence 
donor. 

The blackbod y temperature and radius (scaled to a dis- 
tance of 6.5 kp dReig et aT]|2004l) of SAX J2 103. 5 +454 shown 
in Fig. [9] (open circle) were obtained during a bright state 
dlnam et al.ll2004l) . which may explain the isolated position in 
the plot. Note that the value of the radius and blackbody tem- 
perature are consistent with those of the low-luminosity Be/X- 
ray binaries (filled circles). The difference appears in the X- 
ray luminosity. The higher luminosity in SAX J2103. 5+454 can 
be explained if an accretion disc was present during the obser- 
vatio ns. The detectio n of a quasi-periodic oscillation at 0.044 
Hz (llnam et al.l 2004) an d the corre lation between spin-up/down 
with X-ray flux dBavkal et al.l 12007b indeed supports the pres- 
ence of an accretion disc. 

6. Conclusion 

We have finally unveiled the nature of the compact companion in 
the high-mass X-ray binary 4U 2206+54. The RXTE light curve 
of 4U 2206+54 appears to be modulated with a period of ~1.54 
hours that we interpret as the spin period of the neutron star. 
Evidence for this modulation is also found in INTEGRAL and 
EXOSAT data. Thus we rule out that the periodicity is due to the 
RXTE orbit. 4U 2206+54 becomes the slowest rotating neutron 
star in a high-mass X-ray pulsar with a main-sequence compan- 
ion. The X-ray flux and the spectral index of the power-law com- 
ponents that fit energy and power spectra exhibit variability on 
time scales of days. The amplitude of variability of the X-ray 
flux changed by a factor of 5 on time scales of days. The power- 
law index is larger (i.e. steeper spectrum) during low-flux states 
in the energy spectra and during high-flux states in the power 
spectra. We attribute this variability to changes in the mass ac- 
cretion rate as the neutron star orbits the O-type companion in a 
moderately eccentric orbit. The time span covered by the obser- 
vations is shorter than any of the two suggested orbital periods, 
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